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Introducción 


Cuando el marcapaso natural del corazón no funciona bien y se produce una arritmia, es posible 
restituir el ritmo cardíaco a valores adecuados utilizando los impulsos eléctricos de un marcapa- 
so artificial. 

Como la estimulación eléctrica del corazón presenta un umbral, si se fija una duración para el 
pulso estimulador, su amplitud debe superar un valor mínimo para que el equipo pueda capturar 
el miocardio, provocando un potencial de acción que se propague a toda la masa de músculo y 
pueda desencadenar una contracción ventricular adecuada para mantener la función de la bomba 
cardíaca 

Los marcapasos constan de un generador de impulsos eléctricos, un cable (o cables) en el inte- 
rior de un catéter de un material aislante y un electrodo (o electrodos) en el extremo del cable (o 
cables). 


Una parte de los problemas de investigación que se plantean a propósito de la estimulación eléc- 
trica crónica del corazón involucran directa o indirectamente a los electrodos de marcapaso. En 
el momento actual (comienzos de la década del ochenta) la migración de los electrodos y las 
variaciones temporales del umbral de excitabilidad eléctrica del corazón, son dos de los princi- 
pales problemas que debe afrontar el tratamiento por medio de marcapasos artificiales. 

Para mantener un margen de seguridad razonable es necesario estimular el miocardio con equi- 
pos que liberen, en cada pulso, una cantidad de energía considerablemente mayor que la que en 
principio podría ser suficiente para lograr la estimulación. Esto conduce a que las baterías duren 
menos, a que se altere más fácilmente la superficie de los electrodos y a que el costo del equipo 
sea sensiblemente mayor. 

La aparición de equipos con salida regulable ha venido a aumentar, por su parte, el interés que 
presenta el estudio del así llamado “incremento fisiológico” del umbral de excitabilidad del 


1 El presente trabajo es una versión modificada de “Aspectos Biomatemáticos de la Cardioestimu- 
lación Eléctrica Crónica”, publicado originalmente en abril de 1983. Ese libro describe con cierto 
detalle un trabajo de investigación presentado por el autor, el 18 de agosto de 1982, en el I Congre- 
so Internacional de Biomatemática, realizado en San Miguel de Tucumán, Argentina. Del texto 
original, en la presente versión se suprimieron los detalles del tratamiento electroquímico de la 
interfaz entre los electrodos y los tejidos biológicos. En compensación se introdujo un grado mayor 
de detalle en el tratamiento del umbral y se abordan los umbrales para pulsos de voltaje controlado 
con un grado mayor de generalidad, siguiendo los trabajos “Teoría de Electrodo de Campo Cons- 
tante en la Cardio-estimulación Eléctrica”, presentado en el II Congreso Internacional de Bio- 
matemática, llevado a cabo en Buenos Aires, Argentina, 1984 y “The time constants for cathodic 
make stimulation of electrical syncytia: an application to cardiac pacing”, publicado en Proceedi- 
ngs of the Annual Conference of the IEEE Engineering in Medicine and Biology Society, Philadel- 
phia, USA, 2006. Pese a sus limitaciones (en 2019 bien evidentes) no se actualizó la totalidad del 
contenido debido a que, en opinión del autor, el trabajo original presenta algunos aspectos que aún 
hoy, no son bien conocidos o incluso, al parecer, no han sido publicados por otros. Los avances más 
recientes se mencionan en un Apéndice y los trabajos correspondientes se pueden encontrar en la 
bibliografía complementaria. 
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corazón, desde que un conocimiento lo suficientemente preciso de la evolución temporal del 
umbral, adaptado a cada paciente, permitiría ajustar la salida del equipo estimulador regulable 
hasta un mínimo compatible con la seguridad en el tratamiento. 


En un trabajo previo (que en parte es una base experimental para el presente trabajo) se estudió 
la evolución de las variables del umbral que se utilizan en la práctica clínica para caracterizar el 
umbral de excitabilidad del corazón humano, se analizó el origen de las modificaciones que 
sufre la impedancia del cátodo estimulador y se construyó un predictor estadístico del umbral a 
partir de los valores medidos en el momento de la implantación y durante los primeros tres días 
a partir de ese momento (Fiandra, Suárez Antola y otros, 1981). 

En lo que sigue se comenzará por introducir algunos conceptos (como la distinción entre varia- 
bles locales y variables globales del umbral) previos a la construcción de un modelo matemático 
destinado a contribuir a elucidar algunos puntos delicados de la conexión entre las variables 
medidas en la práctica clínica y los parámetros que caracterizan al sistema electrodo-miocardio. 
Luego se construirá un primer modelo matemático conceptual del sistema electrodo-tejidos no 
excitables interpuestos-miocardio. 


Con el fin de facilitar el trabajo del lector, recordemos algunos aspectos generales de los marca- 
pasos cardíacos artificiales, enfatizando lo referente a los electrodos. 

Cuando el electrodo en el extremo es único, el circuito eléctrico se cierra, a través del conductor 
de volumen de los tejidos, entre ese electrodo que actúa como cátodo y un ánodo formado por 
parte de la superficie externa del generador de pulsos. 

La correspondiente modalidad de estimulación se denomina unipolar. 

Cuando el catéter presenta, además del electrodo en su extremo (cátodo), otro electrodo próxi- 
mo al primero (ánodo), el circuito eléctrico a través de los tejidos se cierra entre esos dos elec- 
trodos. En este caso, la superficie externa del generador de pulsos se encuentra completamente 
aislada desde el punto de vista eléctrico. 

La modalidad de estimulación se denomina bipolar. 

Se puede denominar derivación al sistema formado por el electrodo endocavitario, un electrodo 
indiferente integrado a la superficie de la caja que encierra la batería y los circuitos electrónicos 
de estimulación y detección, y la conexión eléctrica entre ambos electrodos a través de los cir- 
cuitos mencionados. 

La Figura 1 muestra un esquema de un marcapaso con un único electrodo próximo al miocardio 
(electrodo unipolar), ubicado en el interior del ventrículo derecho, apoyado en el endocardio 
(estimulación endocárdica). 
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Figura 1 Esquema de un marcapaso con electrodo unipolar. La cabeza del electrodo actúa como 
cátodo y el generador de pulsos implantado actúa como ánodo. 
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La Figura 2 presenta un esquema del catéter con el cable conductor en su interior, la cabeza del 
electrodo junto con elementos de fijación pasiva al miocardio en un extremo y el conector con el 
generador de pulsos en el otro extremo. 
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Figura 2 Esquema de un catéter con la cabeza de electrodo en uno de sus extremos y un conector 
para el generador de pulsos en el otro. 


El catéter con el cable se puede hacer ingresar a través de la vena subclavia izquierda, entrando 
en la aurícula derecha a través de la vena cava superior, atravesando la correspondiente válvula 
aurículo-ventricular y alojando su extremo en el vértice anfractuoso del ventrículo derecho, 
donde los elementos de fijación del extremo del catéter contribuyen a mantenerlo en posición. 
El par de electrodos (el electrodo unipolar en el extremo del catéter y el electrodo situado en la 
superficie externa del generador de pulsos) se emplea tanto para estimular el corazón como para 
medir la actividad eléctrica espontánea del miocardio. 


En la estimulación crónica el generador de pulsos produce un tren de pulsos de voltaje controla- 
do entre el ánodo y el cátodo, con una separación temporal entre pulsos y una duración de pulso 
adecuadas para el caso. Si se fija una duración para el pulso estimulador de voltaje controlado, 
su amplitud debe superar un valor mínimo (umbral) para que el equipo pueda capturar el mio- 
cardio, provocando un potencial de acción que se propague a toda la masa de músculo. El um- 
bral se suele determinar aplicando el pulso fuera de los períodos refractarios absoluto y relativo 
(umbral diastólico tardío). A menor duración del pulso, mayor resulta la amplitud umbral. 


Si U, representa la amplitud umbral de un pulso rectangular de voltaje controlado y Í pes su 


duración, la dependencia entre ambos se describe a menudo en la práctica clínica mediante la 
fórmula hiperbólica de Weiss-Lapicque: 


U, =U p | 1+ tera (1) 
Pp 
En esta expresión el umbral se caracteriza mediante dos parámetros: la amplitud de voltaje reo- 
básica U z, (la amplitud umbral mínima, por debajo de la cual por larga que sea la duración del 


pulso estimulador, la captura ya no es posible) y la cronaxia /,, ¡; (duración del pulso rectangu- 


lar de voltaje para la cual la amplitud umbral es el doble de la reobase). 

De acuerdo con la fórmula (1), la amplitud umbral disminuye y se aproxima a la reobásica a 
medida que la duración del pulso aumenta tendiendo a infinito. En la práctica la reobase de vol- 
taje se alcanza para una duración finita (duración útil hasta la reobase de voltaje) que a menudo 
está próxima a los 2 ms. A menudo se obtienen valores de décimas de ms, generalmente valores 
comprendidos entre 0.2 y 0.5 ms, para las cronaxias de voltaje. Las reobases de voltaje varían 
más, pero en general toman valores de voltios o fracciones de voltio. 


También se determinan los umbrales de pulsos de corriente controlada para diferentes duracio- 
nes de pulso. En la práctica clínica los resultados de las mediciones del umbral de corriente se 
suelen expresar también recurriendo a una fórmula de Weiss-Lapicque: 


Aspectos Biomatemáticos de la Cardioestimulación Eléctrica Crónica 
(Biomathematical Aspects of Chronic Electrical Cardiostimulation) 
April 1983 (reviewed July 1919) 


lo, 
y O (0 (2) 
Pp 
En esta fórmula /,, es la reobase de los pulsos rectangulares de corriente controlada (que en la 
práctica también se alcanza para una duración finita, la duración útil hasta la reobase de corrien- 


te) y Í,, es la correspondiente cronaxia de corriente. Para un mismo sistema electrodo- 


miocardio, cuando se las determina en el mismo momento, por lo general la cronaxia de voltaje 
resulta algo menor (o a lo sumo igual) que la cronaxia de corriente. 


En trabajos que describen los resultados de las determinaciones de curvas amplitud umbral ver- 
sus duración para pulsos rectangulares de corriente controlada, más orientados hacia la electro- 
fisiología cardíaca, se utiliza a veces la fórmula de Lapicque-Hill. En este último caso el umbral 


de corriente se caracteriza mediante la misma reobase y la denominada constante de tiempo +, ,, 


a través de la expresión: 
1, =—*— (3) 


]=e 


La constante de tiempo 7, , es proporcional a la cronaxia: 


£,, =[(In2)-1 (4) 


cri 


Para una misma duración del pulso, los resultados experimentales aportan evidencia acerca de 
que los umbrales medidos dependen, entre otras cosas de los siguientes factores: 


(1) La ubicación del electrodo en el interior de la cámara cardíaca estimulada. 

(Q) El estado del miocardio excitable, influido por metabolitos, señales endócrinas y neurales, 
globales y locales. Los ritmos circadianos, el sueño, el consumo de alimentos, la actividad físi- 
ca, la fiebre, las drogas, las emociones y otros factores inciden sobre los valores del umbral. 

(3) La reacción inflamatoria desencadenada por la presencia del electrodo, que termina con la 
formación de un casquete de tejido fibroso no excitable en torno a la cabeza del electrodo. Esta 
reacción es la principal responsable del “incremento fisiológico” del umbral. ? 


(4) La mayor o menor rigidez del catéter. * 


(5) Las propiedades del electrodo que se encuentra en su extremo (tamaño, forma, material y 


2 En los capítulos de Doenecke y coautores (Doenecke, Flóthner y otros, 1975) y de Szabo y Solti 
(Szabo y Solti, 1975) se pueden hallar resultados sobre la evolución de las variables del umbral y en 
particular sobre el incremento fisiológico. 


3 El grado de rigidez del catéter tiene influencia sobre la reacción inflamatoria desencadenada por 
el electrodo y sobre las fluctuaciones en las variables del umbral medidas en cada paciente en un 
instante dado, es decir, para un mismo estado del miocardio. Esas fluctuaciones, que suelen ser 
significativas cuando se emplea un electrodo colocado en el extremo de un catéter rígido (como los 
destinados a estimulación transitoria) disminuyen mucho cuando el electrodo se encuentra ubicado 
en el extremo de un catéter complaciente (como los destinados a la estimulación crónica). 
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textura). * 
(6) La distancia entre el electrodo y el tejido excitable. 
(7) El procedimiento utilizado para determinar el umbral. ? 


En consecuencia, las curvas amplitud umbral versus duración para pulsos de corriente y de vol- 
taje controlado son una propiedad del sistema electrodo-miocardio considerado como un todo. 
Describen lo que puede denominarse excitabilidad global de ese sistema. 


La excitabilidad global debe distinguirse de la excitabilidad local. Esta última es la excitabili- 
dad definida a nivel de las membranas de las fibras miocárdicas donde se produce el patrón de 
polarización originado por el pulso estimulador. 


Para discutir (teniendo en cuenta la distinción entre excitabilidad global y excitabilidad local) la 
influencia del tamaño, la forma, el material y la textura de los electrodos de marcapaso sobre los 
umbrales de captura del miocardio, es conveniente comenzar por una descripción funcional del 
sistema electrodo-miocardio. 


Descripción funcional del sistema electrodo-miocardio. 
Los electrodos y los tejidos biológicos del paciente forman un sistema cuyas propiedades eléc- 
tricas, electroquímicas, mecánicas y biológicas se deben tener para estudiar la relación entre las 


variables globales y las variables locales del umbral. 


Si se enfatiza el flujo de señales a través del sistema *, se pueden distinguir cuatro subsistemas 
parcialmente superpuestos en el espacio. Ver Figura 3. 


4 Durante la medición de las variables del umbral efectuadas con el electrodo de auto-fijación in- 
ventado por Orestes Fiandra (Fiandra, Espasandín y otros, 1976) se halló evidencia acerca de la 
influencia de la forma del electrodo sobre el umbral. La influencia del tamaño del electrodo y del 
espesor del casquete de tejidos no excitables interpuesto entre la superficie del electrodo y el mio- 
cardio excitable también es bien conocida. Cuando se mide el umbral con pulsos de voltaje contro- 
lado, es conocida la influencia de la polarización de la interfaz electro-tejidos, cuya importancia 
sobre el umbral puede depender del material del electrodo y de su textura (Irnich, 1975; Irnich, 
1980). 


5 Cuando el umbral se mide comenzando por pulsos supra-umbrales cuya amplitud se disminuye 
hasta que el marcapaso pierde la captura del miocardio, los umbrales así determinados resultan 
inferiores respecto de los valores que se obtienen cuando el umbral se mide comenzando por pulsos 
sub-umbrales cuya amplitud se va incrementando hasta que el equipo logra capturar la actividad 
eléctrica del miocardio (Christer, Sylvén, y otros, 1982). 


$ Los fundamentos del análisis de sistemas fisiológicos basado en el flujo de señales pueden hallarse 
en el libro de Talbot y Gessner (Talbot y Gessner, 1973). 
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Figura 3: Esquema para una descripción funcional del sistema electrodo-miocardio cuando se enfa- 
tiza la estimulación eléctrica del corazón. 


(A)-Un subsistema eléctrico pasivo, que denominamos área bioeléctrica pasiva. Está compues- 

ta por cuatro módulos conectados en serie: 

(1)El catéter con el cable y el electrodo activo (cátodo) en su extremo proximal al miocardio (a 
lo que se añade el cable del ánodo y el ánodo si la estimulación es bipolar). 

(Q)La interfaz entre la superficie del electrodo activo y los tejidos biológicos. 

(3) Un casquete de tejidos no excitables que se interpone entre la superficie del cátodo y el mio- 
cardio excitable, en la fase crónica de la estimulación consiste en una cápsula de tejido fibro- 
so. 

(4) El conductor de volumen formado por el resto de los tejidos, hasta el electrodo indiferente. 


En el área bioeléctrica pasiva se pueden distinguir dos señales de entrada y una señal de salida. 
Una de las señales de entrada es eléctrica: es un pulso de corriente que atraviesa cada uno de los 
cuatro módulos situados entre el electrodo activo y el electrodo indiferente. La otra señal de 
entrada es cinemática, proviene del cuarto subsistema y será considerada luego de describir 
dicho subsistema. La señal de salida es una distribución espacial de corriente eléctrica que co- 
necta este subsistema con el subsistema siguiente. 


(B)-Un subsistema eléctrico activo, que denominamos área bioeléctrica activa. Está compuesto 
por tres módulos, con los dos primeros conectados en serie: 

(1)El tejido cardíaco excitable situado próximo al cátodo, blanco de la estimulación. 

(2) El resto del miocardio. 

(3)La jerarquía de marcapasos naturales. 


La señal eléctrica de salida desde el área bioeléctrica pasiva (el patrón de distribución de la den- 
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sidad de corriente eléctrica en el conductor de volumen de los tejidos) es la señal de entrada al 
área bioeléctrica activa. En el área bioeléctrica activa se pueden distinguir dos señales de salida: 
un patrón de ECG cuyas características permiten constatar si la captura del miocardio se ha 
producido o no, y una señal fisicoquímica (potencial de acción) que conecta este subsistema con 
el siguiente. 


(C)-Un subsistema mecánico activo o área biomecánica activa. Formado por el músculo car- 
díaco contráctil, recibe la señal fisicoquímica del área bioeléctrica activa y se contrae produ- 
ciendo una señal mecánica consistente en un campo variable de esfuerzos que lo conecta con el 
subsistema siguiente. 


(D)-Un subsistema mecánico pasivo o área biomecánica pasiva. Incluye los tejidos que rodean 
la cabeza del electrodo activo y el catéter, considerados desde el punto de vista mecánico. Reci- 
be la señal mecánica proveniente del área biomecánica activa y produce una señal de salida 
cinemática. 

Esta señal cinemática entra en el área bioeléctrica pasiva donde se evidencia en una modifica- 
ción del espesor de la cápsula de tejido no excitable interpuesto entre la superficie del electrodo 
activo y el miocardio excitable blanco de la estimulación. 

Se cierra así un circuito de transferencia de señales, que comenzando con la entrada del pulso de 
corriente en el área bioeléctrica pasiva termina con una señal cinemática de entrada a esta mis- 
ma área, proveniente del área biomecánica pasiva. 


Este bucle de retroalimentación complica la construcción de un modelo matemático del sistema 
electrodo-miocardio durante la estimulación por electrodos de marcapaso, debido a que la dis- 
tancia entre el electrodo y el miocardio excitable puede variar de un pulso a otro, dependiendo 
del grado de complacencia del catéter, de la frecuencia del tren de pulsos estimuladores, de la 
amplitud de la contracción del músculo cardíaco, de la ubicación y anclaje del electrodo, entre 
otros parámetros. 

No obstante, es posible abrir el lazo de retroalimentación considerando la distancia electrodo- 
miocardio excitable como una variable aleatoria, que toma un valor constante durante la dura- 
ción de un mismo pulso estimulador, aunque ese valor pueda modificarse de un pulso a otro. 
Esta constancia para un mismo pulso se debe a que los pulsos poseen una duración del orden de 
las milésimas de segundo, mientras que las fases del ciclo cardíaco y los fenómenos mecánicos 
en el sistema electrodo-miocardio se desarrollan en escalas del orden de las décimas de segun- 
do. 


Una vez que se implanta el electrodo, la presencia del material extraño produce una reacción 
inflamatoria seguida de la formación de una cápsula de tejido fibroso que rodea la cabeza del 
electrodo. Como consecuencia, las propiedades locales del conductor de volumen cambian y 
debido a este cambio cabe esperar que puedan cambiar las propiedades de detección de señales 
y de estimulación del electrodo considerado como parte integrante del sistema electrodo-tejidos 
no excitables interpuestos-miocardio. 


Un modelo de cable tridimensional excitable 


En su libro sobre el flujo de corriente en los tejidos excitables, Jack, Noble y Tsien (1975) des- 
criben las membranas excitables de un sincicio eléctrico (como el miocardio) empleando un 
cable no lineal homogéneo e isótropo, que llena el espacio, de tal forma que en cada punto se 
superponen el espacio extracelular, el espacio intracelular y la membrana excitable. 

Se sabe que el conductor de volumen formado por el miocardio no es isótropo y tampoco es 
homogéneo (Trautwein, 1975; Chapman y Fry, 1978). Además, hay evidencia acerca de que el 


Y 
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grado de anisotropía del conductor de volumen intracelular es diferente al grado de anisotropía 
del conductor de volumen extracelular (Chapman y Fry, 1978). 

Las uniones entre células del miocardio ventricular presentan una resistencia baja al flujo de 
corriente eléctrica (Figura 4) por lo cual se justifica postular un continuo intracelular. 
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Figura 4: A la izquierda un esquema de las uniones entre células miocárdicas a través de discos 


intercalares. Estos discos presentan nexos atravesados por canales iónicos gigantes que comunican 
entre sí a las células adyacentes y contribuyen a la unidad de ambiente iónico en el espacio intrace- 
lular. A la derecha un esquema de la disposición de las fibras miocárdicas en relación con los capi- 
lares. 


Si representamos por Í el instante de tiempo considerado y por Y el vector de posición del 
punto del espacio considerado, asumiendo homogeneidad pero no isotropía, es posible generali- 
zar la ecuación de Jack, Noble y Tsien, postulando la siguiente relación entre la densidad de 


corriente transmembrana total /,, (1 $ ) por unidad de volumen (corriente que se puede expre- 


sar como suma de una corriente de conducción a través de canales iónicos 1, lv, w) y de 
una corriente de desplazamiento C,,, EM relacionada con la variación en la carga eléctrica 


en el capacitor de membrana), el campo de voltaje transmembrana V,, (1,7 ) y la densidad de 


corriente exógena /., (E F) , por unidad de volumen: 


» 1 0 lo 0? A 
Eo MA a A 6) 
"Ot R, Ox R, 0y R. Oz 


Xx 


El voltaje transmembrana V,, se define como la diferencia entre el potencial eléctrico en el es- 


pacio intracelular (formado por las células eléctricamente interconectadas) menos el potencial 
eléctrico en el espacio extracelular (intersticial). 
El parámetro C, y es una capacidad eléctrica de membrana por unidad de volumen. 


El vector W representa las variables de activación, recuperación y adaptación asociadas a los 
canales iónicos de la membrana excitable. 


Los parámetros RR, ,R. representan ciertas resistencias eléctricas por unidad de longitud en 


dirección de los ejes X,y,Zz de un sistema de coordenadas cartesianas ortogonales. 


En un conductor lineal y anisótropo, la densidad de corriente y el campo eléctrico se 
relacionan a través de un tensor de conductividad simétrico (Landau y Lifchitz, 1960). 
Aquí suponemos que las direcciones de los ejes de coordenadas coinciden con las direcciones 
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principales del tensor. Por ese motivo son suficientes tres parámetros para describirlo. 


Si bien la ecuación (5) representa un cable tridimensional anisótropo, no tiene en cuenta la an- 
isotropía desigual que presenta el miocardio. 

Tampoco tiene en cuenta que las direcciones principales del tensor de conductividad deben ro- 
tan a medida que se consideran capas de músculo ventricular cada vez más profundas respecto 
del endocardio. Cabe esperar este comportamiento del tensor de conductividad debido a que la 
dirección de las fibras del miocardio ventricular rota al pasar de una capa a otra (Burton, 1972; 
Davson y Eggleton, 1968; Urschel y Sonnenblick, 1974). 


Las corrientes eléctricas se toman positivas cuando fluyen desde el espacio intracelular hacia el 
espacio extracelular. 

Las derivadas parciales que aparecen en la ecuación (5) se toman en la dirección de los 
ejes coordenados. 


Suponemos un estado cuasi estacionario tal que la densidad de corriente exógena, por 
unidad de volumen, se puede factorizar así: 


1,Lt,F)=5(0): EJ(5) (6) 
En la formula (6) ¡(e ) representa la corriente eléctrica inyectada a través del electrodo próximo 


al miocardio, mientras que F, (7) es un campo escalar que describe la distribución espacial de 


la corriente inyectada en el conductor de volumen de los tejidos. 


Si la duración del pulso estimulador es de un orden de magnitud mayor que la escala de tiempo 
de los canales iónicos asociados con las variables de activación de la membrana, podemos su- 
poner que en cada instante de tiempo las variables de activación se encuentran relajadas al equi- 
librio con el voltaje transmembrana. 

Asimismo, si la duración del pulso estimulador es de un orden de magnitud inferior al de la 
escala de tiempo de las variables de recuperación, podemos suponer que estas últimas se en- 
cuentran congeladas en sus valores correspondientes a la membrana en reposo. 

En el caso del miocardio, estas hipótesis simplificadoras pueden ser aceptables, al menos en una 
primera aproximación (Fozzard, 1979). 


Entonces la densidad de corriente iónica /, (v, ) queda en función del voltaje transmembrana 


lon 


y se comporta como sugiere la Figura 5. 


Figura 5: Gráfica cualitativa de la corriente iónica transmembrana cuando las variables 
de activación están en equilibrio con el voltaje transmembrana y las variables de recupe- 
ración se encuentran congeladas en sus valores de reposo. 
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La corriente presenta tres ceros: el situado a la izquierda es el voltaje de reposo V,, el 
intermedio es el denominado umbral uniforme de membrana V,, y el último cero V, , situado a 


la derecha, aparece como consecuencia de haber congelado las variables de recuperación (Foz- 
zard, 1979; Jack, Noble y Tsien, 1975). 

Teniendo todo esto en cuenta, la ecuación (5) del cable anisótropo, tridimensional y no lineal se 
puede escribir así: 


Pe E A E 
m m 
RON R, dy R. 02 
Suponiendo una polarización espacialmente uniforme de las membranas del cable tridimensio- 
nal, la ecuación a derivadas parciales (7) se reduce a una ecuación diferencial ordinaria. Una vez 
que la corriente inyectada a través del circuito del cátodo se anula, la dinámica del voltaje 
transmembrana viene dada por la ecuación no lineal: 


d 
Ca ] e 5 Ea ja (8) 
y Ev, == (V,) 
De los tres ceros de /, (v,. ), Vi y V,son puntos de estado estacionario estable de esta última 


lon 


V, +i(t)-F(F) (0) 


ecuación, mientras que V,, es inestable. 


Si la polarización uniforme lleva a V,, (1 ) desde V, hasta un valor inferior a V,, , el voltaje 


m 


Vv (1) retorna directamente a su valor de reposo V, . 
Si la polarización uniforme lleva a V,, (1 ) desde V hasta un valor superior a V,, , se produce la 


excitación (V,, (1 ) tiende a V,, ). Por tanto, V,, es efectivamente el denominado umbral uniforme 


de membrana. 
S1 la polarización producida no es espacialmente uniforme, una vez cesado el pulso de corriente 


estimuladora, el voltaje transmembrana V (r, r ) evoluciona a partir de una distribución inicial 


m 


Y, (0, r ) de acuerdo con la ecuación a derivadas parciales no lineal siguiente: 


PA A A Ai A (9) 


Elipsoides críticos de células miocárdicas y sus dimensiones. Regiones de influencia 
del electrodo y familias de estados de polarización umbral. 


En el caso de una fibra de Purkinje, para producir un potencial de acción mediante una 
polarización no uniforme provocada por un microelectrodo actuando como cátodo puntual, una 
longitud mínima de membrana de la fibra debe ser despolarizada por encima del umbral 
uniforme de membrana (Fozzard and Schoemberg, 1972). 

Esta longitud mínima corresponde a lo que Rushton denominó “liminal length” en el caso de 
una fibre nerviosa (Jack, Noble y Tsien, 1975). 

Noble (Jack, Noble y Tsien, 1975) dedujo la siguiente fórmula para la longitud límite de una 
fibra excitable estimulada mediante un microelectrodo (cátodo puntual): 


Al, (Vi, ) 


10on 


dv, 


m 


lim 
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1Oon 


En esta expresión =—— = 
R dv, 


m m 


- (0, V,) b(v, VA Y +c(V,,—V,)' | siendo b y C paráme- 


m 


d IL, (V,) 


y 4,, es la constante de espacio de la fibra. 


Noble asumió /,,, (v,, ) = 
tros positivos. Obtuvo un único perfil umbral, estacionario, para el voltaje transmembrana en 
un modelo de fibra de extensión infinita. 


El miocardio ventricular sobre el cual actúa el cátodo no se comporta como una fibra, sino como 
un sincicio eléctrico. Por ese motivo se plantea aquí un modelo basado en un cable 
tridimensional. 

Los resultados experimentales presentados por Lindemans en su tesis de doctorado (Lindemans, 
1977) sugieren que cuando el miocardio ventricular es estimulado por un cátodo externo, es 
necesario despolarizar una masa crítica de células excitables por encima del umbral uniforme de 
membrana para que pueda emerger un potencial de acción capaz de propagarse y desencadenar 
la contracción del músculo ventricular. 


Partiendo de la ecuación (9), un análisis heurístico permite obtener una expresión aproximada 
para las dimensiones de esa masa crítica de células. 
Se puede suponer que el electrodo produce patrones elipsoidales de puntos con igual 


polarización. Representemos por ¡0 y and l z a las distancias entre el punto más despolari- 


zado del miocardio ventricular (despolarizado por encima del umbral uniforme de membrana) y 
los puntos ubicados sobre los ejes de coordenadas Xx, Y ,Z donde el voltaje transmembrana es 


igual al umbral uniforme de membrana V,,, . 


Estas distancias caracterizan una región de influencia del cátodo sobre el tejido excitable. 


En el punto donde la membrana está más despolarizada aproximamos así a las derivadas parcia- 
2 

o” (v,, Vi) O Voz (v, Vir) o? (v, Vr) 

A A En 2 

Ox l 19) y l A OZ l 

El fundamento de este tipo de aproximación se puede hallar, por ejemplo, en el libro de Kyrala 

(Kyrala, 1967). 


Q 


les: 


En ese mismo punto aproximamos:  Í,,, (v, ) = (E Í, (Ve ) (v, — Ve, ) 


0 : : do 
Teniendo en cuenta que na Le (Vr, ) es negativa, sustituyendo todas estas aproximaciones en 
m 


la ecuación (9) obtenemos, finalmente, para el punto más despolarizado: 


O 
O v' UN = 
"Ot 
1 1 
+= += V, 
RECREO RA nm) 
polarizado tiende a disminuir y este resultado sugiere que el miocardio estaría en un estado 
subumbral. 


a hn) : : E Jo.) (1D 


> + + 
ov, RL R, E RL 


Si 


es mayor que el voltaje en el punto más des- 


10) 
OV A 


m 
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1 1 1 


Si + —— + 
O A A 


es menor que |—1,, VA ) el voltaje en el punto más despo- 
z z oV,, 

larizado tiende a aumentar y ahora el resultado sugiere que el miocardio estaría en un estado 

supraumbral. 

Finalmente, si son iguales, ese voltaje se mantendría constante en el nivel de aproximación en el 

cual nos encontramos y esto sugiere que el miocardio se encontraría en un estado umbral. 


Admitiendo estas suposiciones, resulta que una región crítica puede ser considerada como un 
mi Ec e A a 
elipsoide de semidiámetros 1,” , 1,” y Ll” que verifican: 


1 1 1. al, (Vo) 


(12) 


De (12) se desprende que no hay un único elipsoide crítico, puesto que solo establece una 
restricción sobre los valores de los semidiámetros, pero no los determina, a diferencia de lo que 
ocurre en el caso de la fibra, donde la longitud límite queda determinada por los parámetros 
electrofisiológicos y geométricos.” 

Por tanto no hay un único perfil de polarización umbral en el caso del cable tridimensional no 
lineal, sino toda una familia de estados umbral. 


La excitabilidad global (determinada mediante las variables del umbral que se utilizan en 
clínica) se puede relacionar con la excitabilidad local (determinada por los umbrales de 
membrana medidos en experimentos electrofisiológicos) a través de esas regiones críticas de 
células miocárdicas, cuyas dimensiones se pueden estimar a partir de la ecuación (12). 


En principio, de las ecuaciones que describen el flujo de corriente eléctrica en el conductor de 
volumen de los tejidos y generan la polarización de las membranas excitables, acopladas con las 
ecuaciones del cable tridimensional no lineal que describen la dinámica de las membranas, se 
podría obtener una condición que deben cumplir los pulsos de corriente eléctrica (e) 


introducidos en el área bioeléctrica pasiva para poder capturar la actividad eléctrica del 
miocardio justo en el umbral. 


Si esta condición se aplica a un tren de pulsos lo bastante separados como para que la influen- 
cia del pulso precedente no sea significativa, es suficiente considerar un único pulso aislado 
para analizar los fenómenos relacionados con el umbral. 

Por lo general, éste es el caso en la estimulación mediante electrodos de marcapaso. 

El pulso de corriente puede ser tanto un pulso de corriente controlada como el pulso de corriente 
asociado con un pulso de voltaje controlado. 


Deducción preliminar de la fórmula de Lapicque-Hill: análisis de los umbrales de 
captura mediante pulsos de corriente controlada. 


7 Si en la fórmula (12) se hacen tender / * y Le a infinito, se multiplican ambos miembros por la 


. . . R, 
resistencia de la membrana en reposo R,, , y se tiene en cuenta que R = AL es una constante 


Xx 
de espacio en dirección del eje X, despejando la dimensión crítica de , se obtiene la misma fórmula 


hallada por Noble (fórmula (10)) para estimar la longitud límite de una fibra excitable. 
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Aún en el caso en el cual se asuma que las variables de activación se encuentran relajadas al 
equilibrio con el voltaje transmembrana y las variables de recuperación se encuentran 
congeladas en sus valores de reposo, el problema matemático es, al parecer, lo bastante 
formidable como para, posiblemente, requerir un abordaje numérico ab-initio. 

No obstante, añadiendo un conjunto de hipótesis simplificadoras muy restrictivas, es posible 
intentar un abordaje analítico y deducir la fórmula de Lapicque-Hill para pulsos rectangulares 
de corriente controlada, o sea la ecuación (3), que re-escribimos aquí para mayor comodidad del 
lector: 


1 
A G bis) 


Comenzando nuevamente con la ecuación (9), y siguiendo la línea de razonamiento esbozada en 
el Apéndice, se puede obtener la siguiente condición para que un pulso de corriente catódica 


¡(e ) de duración 1, sea supraumbral y genere un potencial de acción: 


o Ñ 
[76, (1, -1)i(e)ar>0, (13) 
Un pulso es justo umbral si la desigualdad (13) se convierte en una igualdad y es subumbral si la 
desigualdad se invierte. En esa desigualdad (Y. es una carga umbral. 


La función sin dimensiones G, (1) es una exponencial decreciente, definida para 1> 0: 
Í 
aloe E (14) 
s,Í 


Si el pulso de duración 1, es de amplitud constante /,, , la condición para que el pulso sea justo 


umbral se reduce a: 


1 


t, , E 
de G, (1, -1)i(t)dt=1,1, | 1-e ** |=0, (15) 
De (15) se deduce (3 bis) con la reobase de corriente: 
O; 
Exp = a (16) 


s,Í 
La constante de tiempo en la fórmula de Lapicque-Hill que relaciona la amplitud umbral de un 
pulso catódico con su duración, viene dada a su vez por la fórmula: 
2 
(1 +e) 


ey AE (17) 
ds yA +(n+e) 


En la fórmula (17) los parámetros que allí aparecen tienen la siguiente interpretación: 


- 1, es el radio del electrodo, que puede definirse así, en función del área A, de la cabeza del 


electrodo: 1, = A 
T 


- € es el espesor del casquete de tejidos no excitables interpuesto entre la superficie del electro- 
do y el miocardio excitable. 
Ñ 1% = Ro, Co 


y es la constante de tiempo de la membrana. 
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Ls dub 1 R R 
- 2 =5 | 7 +57 | siendo, por definición A. == y ds =5 
A NA A R, R, 


- y es un número positivo sin dimensiones (del orden de la unidad) que resulta de aproximar la 
2 1 2 


. + . 
Ry 0 7 RON 


Xx 


expresión y en los planos que cortan ortogonalmente al eje z,em- 


pleando la técnica de aproximación ya mencionada (Kyrala, 1967): 


2 2 2 
En Í z 2 Y + En j z 2 py = Y Ley 2 (v, V,) 48) 
R e Ox R y O y (r, + e) 
En la Figura 6 el miocardio se describe, en coordenadas cartesianas ortogonales, ocupando la 
región 220. 


A — 


¿ER y AA 
ATA co todo O 


| miocardio 
( 
( 


Figura 6. Esquema de un sistema electrodo-miocardio, con un cátodo ubicado a una cierta distancia 
e de la región ocupada por el músculo ventricular. El origen de coordenadas se encuentra en el 
plano que representa la superficie del miocardio. El eje Z es el eje de simetría del cátodo. 


En los puntos situados sobre el eje z la despolarización de las membranas es máxima y decrece 
en cada plano ortogonal al eje z a medida que el punto de la membrana se aleja del eje, como 


4 S gi : X y 
sugiere la figura. Las coordenadas son números sin dimensiones X =-—— Y =-— que se 


xx yy 


introducen porque, también como sugiere la figura, se puede suponer que 
Y, (z, X,Y, z) =V, 6 Vx?+Y?, 2) A medida que VX?+Y? crece para z constante, V,, 


se aproxima al voltaje V, de la membrana en reposo. 
El parámetro Y aumenta cuando las líneas de corriente producidas en el circuito del cátodo se 


separan, como ocurre, en igualdad de las demás condiciones, al pasar de la estimulación mono- 


14 


Aspectos Biomatemáticos de la Cardioestimulación Eléctrica Crónica 
(Biomathematical Aspects of Chronic Electrical Cardiostimulation) 
April 1983 (reviewed July 1919) 


polar a la bipolar. 


J _— ) 
P o ' 
a an od o N | 

p dl 


Y 
DA c3todo 
Fer 
: q + pa 
A AA 
, mi *R ” 7 — > 
sel HA nod 
ph 9 AA m 10 
7 T mircardís / X tar 
- | x A s 
e A 
II? 
: ó / 
Estimv facie» ' Estilmudacion 
bipolar ¿Age ¿bren montpodor ; el ¿nedo 
/ to” ' 
dal Lduras ele corriente Se puede congroderar 


S.tvido e. ef infinito. 


Figura 7. Ilustra la mayor apertura de las líneas de corriente al pasar de la estimulación monopolar 
a la bipolar. 


La fórmula (17) predice que la constante de tiempo f, , debería crecer cuando aumentan el ra- 


dio del electrodo o el espesor del casquete de tejidos no excitables, y debería disminuir, en 
igualdad de las demás condiciones, al pasar de la estimulación monopolar a la bipolar. 

La Figura 8 muestra un esquema de la relación entre la constante de tiempo y la suma del radio 
del electrodo con el espesor del casquete de tejidos no excitables. 


La constante de 
tiempo Es en 
función de r+e 


- - - —> PA 


E*) Fe +e 


Figura 8. La constante de tiempo 7, , en función de 7, +€.. 


La reobase de corriente resulta proporcional a la carga umbral O. y al término 


2 : 
Y A as (r, + e) siendo K una constante de proporcionalidad: 


Ip =K0.[ 142, +(1m+e) | (19) 


Entonces, de la fórmula (19) se desprende que, en igualdad de las demás condiciones, la reobase 
debería aumentar con el radio del electrodo, con el espesor del casquete de tejidos no excitables 
interpuesto entre el electrodo y el miocardio y con el cambio en la modalidad de estimulación de 
monopolar a bipolar. 


A diferencia de otras expresiones que relacionan la reobase con I;, +€ , la reobase dada por la 
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ecuación (19) no tiende a cero cuando 7, +€ tiende a cero. El comportamiento predicho por la 


fórmula (19) parece acorde a los resultados experimentales de Lindemans (Lindemans, 1977). 


En general, para los electrodos convexos convencionales, la evidencia experimental confirma 
estas predicciones cualitativas, tanto en el caso de la constante de tiempo como en el caso de la 
reobase (Irnich, 1975; Irnich, 1980; Doenecke, Flóthner y otros, 1975; Lindemans, 1977). 


Análisis de los umbrales de captura mediante pulsos de voltaje controlado. 


Los generadores de pulsos de los marcapasos utilizan pulsos de voltaje controlado de amplitu- 
des aproximadamente constante (pulsos rectangulares de voltaje). 

Mientras la amplitud del pulso de voltaje entre el cátodo y el ánodo permanece constante, el 
pulso conjugado de corriente eléctrica sufre variaciones que reflejan los cambios en la impedan- 
cia de carga vista por el generador de pulsos. 

Ver la Figura 9, donde se ve la respuesta de voltaje entre el cátodo y el ánodo cuando se inyecta 
un pulso rectangular de corriente controlada, y la respuesta de corriente cuando se aplica un 
pulso rectangular de voltaje controlado entre el cátodo y el ánodo. 


i 


Polarization_—-7 
potential 


Polarization 
potential 


(a) 


Polarization 


Polarization 


(b) 


Figura 9. Formas de pulso conjugadas en el área bioeléctrica pasiva. En (a) se representan dos pul- 
sos rectangulares de corriente controlada de amplitud ¿, junto con los correspondientes pulsos 
conjugados de voltaje U (1) . En (b) se representan dos pulsos rectangulares de voltaje controlado 
de amplitud U, junto con sus pulsos conjugados de corriente i(e). En ambos casos los pulsos con- 


jugados presentan una respuesta óhmica seguida de una respuesta diferida asociada a la polariza- 
ción de los electrodos de marcapaso. 


La impedancia de carga vista por el generador de pulsos es independiente de si el miocardio es 


capturado o no por el marcapaso: refleja las propiedades del área bioeléctrica pasiva solamente. 
Se puede expresar como suma de dos términos. 
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Uno de ellos es una resistencia R p que representa el comportamiento óhmico de la combina- 
ción en serie de los cables del (o los) electrodo (electrodos) representados por una resistencia 
R.. con el conductor de volumen de los tejidos representado por una resistencia Ry denomina- 
da resistencia de los tejidos (Irnich, 1975; Fiandra, Suárez-Antola y otros, 1981). 

La resistencia de los tejidos Ry depende del tamaño y la forma del electrodo estimulador, así 


como de la conductividad eléctrica en su vecindad: no es una propiedad de los tejidos en sí 
mismos (Irnich, 1975; Fiandra, Suárez-Antola y otros, 1981). 


Si U, es la amplitud de un pulso de voltaje aplicado (no necesariamente umbral) e ie) es el 


pulso de corriente conjugado, la resistencia KR, =R +R, se puede determinar así: 


R =lim_.,>2 Por otra parte, si i ¿es la amplitud de un pulso de corriente introducido a 


p 150 i(+) 
E Uli 
través de los electrodos y si U' ( ) es el pulso de voltaje conjugado, R $ lim,_o ut) 
lo 


U, UL) 


En la práctica se suele determinar los cocientes 0) o —— alos 90 ys de iniciado el pulso, 
Í lo 


es decir, para t=90us. 
Debido a que i(e) disminuye más rápido de lo que U (2) aumenta, el valor de R p estimado 


para pulsos rectangulares de voltaje controlado resulta un poco mayor que el valor de R p esti- 
mado para pulsos rectangulares de corriente controlada [24]. Como la resistencia del cable se 
puede hallar en forma independiente, una vez conocida R p Se puede determinar la resistencia de 


los tejidos. En una primera aproximación esa resistencia se puede considerar inversamente pro- 
porcional al tamaño del electrodo. 


El otro término de la impedancia vista por el generador de pulsos es una impedancia que corres- 
ponde a la polarización de la capa doble electroquímica que se forma en la interfaz entre los 
electrodos y los tejidos biológicos (Greatbach y Kahn, 1974; Szabo y Solti, 1975; Irnich, 1975). 
Desde el punto de vista electroquímico las interfaces entre los electrodos de marcapaso y los 
tejidos (o las soluciones fisiológicas) exhiben un comportamiento no lineal. 

Es necesario superar una diferencia de potencial crítica a través de la interfaz para que comience 
a fluir una corriente farádica en paralelo con la corriente capacitiva. 

Una vez superado este valor crítico, la relación entre la corriente farádica y la diferencia de po- 
tencial que la impulsa no es lineal: se parece a lo que se observa conectando en serie dos diodos 
Zener, oponiéndose el uno al otro, como sugiere la Figura 10 para un electrodo de Platino. 


Figura 10. Equivalente eléctrico de la interfaz para un electrodo de Platino. 


Pero para los fines del presente trabajo se utilizará, en una primera aproximación, un modelo 
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lineal análogo al tipo de modelo utilizado para relacionar el esfuerzo mecánico con la deforma- 
ción en la teoría de los materiales hereditarios (Rabotnov, 1980). 


En ese caso las diferencias de potencial U (: ) entre el ánodo y el cátodo se pueden escribir co- 


mo suma de caídas óhmicas en los cables y los tejidos más las caídas de potencial en las capas 
dobles electroquímicas, como se muestra a continuación: 


U(1)=R, 10) [X(=a) ¡(u) du 202) 


Aquí la función de respuesta al impulso K (2) es positiva y caracteriza el comportamiento eléc- 


trico de la interfaz en una primera aproximación lineal. Si R, representa la resistencia de la 


interfaz a la corriente directa: 
| k(+) dr=R, (20 b) 
0 


Las caídas de potencial a través de las capas dobles se toman referidas a la fem de equilibrio 
(Moore, 1962), de modo que esta última se elimina de las ecuaciones. 


Como el generador de impulsos impone U (2 ) , conviene invertir la ecuación (20 a). Se obtiene: 


(9-2 (06) ju JU) de 21a) 


La función L(+) también es positiva y verifica: 


fz(») eE EEN Q1b) 
A R, +R, 


Cuando se aplica un pulso de voltaje de amplitud constante U,, la ecuación (21 a) se reduce a 


la siguiente: 


s O, 
= —— 11 
i(+) 2. +R, (1+p(1) (22) 
Por definición: 
R ROA 
ES L(u) du (23) 
R, R, | 


De (23) y (1 b) se desprende que la función pls) es positiva y monótona decreciente, desde 


R 
p(0) =—2 hasta cero cuando £ —>00. Entonces de (Q2) resulta que el pulso de corriente ile) 
Pp 
disminuye desde el comienzo hasta el final del pulso de voltaje, tal como sugiere la figura 4 (b). 


De (21 c) se desprende que Ll ) tiende a cero más rápido que K (z ) a medida que pasa el 


tiempo. Como consecuencia, en la Figura 9 (b) i (1) disminuye más rápido que lo que aumenta 


U(t) en la Figura 9 (a). 


Un pulso rectangular de voltaje será un pulso justo umbral U, (1) si el pulso de corriente con- 


jugado de duración Í p cumple la condición de umbral establecida empleando el criterio de la 


funcional de excitación para un pulso catódico: 
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t, 


Ja(r, -1) i,(t) dt=0, (24) 
0 
Por otra parte, si U,es la amplitud umbral de un pulso de voltaje de duración Í p» Je la fórmula 


(22) se desprende que: 


, U 
A A (1+p(1) (25) 
Pp d 


De (24) y (5) resulta la siguiente fórmula para la relación entre la amplitud de voltaje y la du- 
ración de los pulsos catódicos: 


a). E 25 


fa(s) ars [6(r,-1) p(t) dt 


t, (1,-:) 
y que ln. fe 3 p(t) dt=0, susti- 
0 


O, 


Teniendo en cuenta: la fórmula (14), que Lo = 

S,I 
tuyendo (14) en (26) se obtiene la siguiente fórmula que relaciona la amplitud umbral de un 
pulso rectangular de voltaje controlado con su duración: 


U 28 
Ú, (r,) = 7 j a (1,1) ( ) 
l-e * l+—|e * p(t) dt 
] a 
En esta última fórmula: 
Uy, =lim, ,. U, (1,)=(R,+R,) 17, (29) 


Así pues, las reobases para pulsos rectangulares de voltaje y las reobases para pulsos rectangula- 
res de corriente se conectan entre sí a través de la resistencia total R mas R, vista por el genera- 


dor de pulsos en condiciones de corriente continua. 

El factor reobase de corriente en la fórmula (29) contribuye a incrementar la reobase de voltaje 
con el aumento en el radio del electrodo, en el espesor del casquete de tejidos no excitables in- 
terpuesto entre el electrodo y el miocardio y con el cambio en la modalidad de estimulación de 


z , 2 2 
monopolar a bipolar, puesto que la reobase es proporcional a y A; ab e + e) z 


No obstante, en el otro factor presente en la mencionada fórmula, aparece R p - La componente 


dominante de este resistor óhmico se comporta en función del radio r, aproximadamente como 


de . Puesto que un factor aumenta y el otro disminuye al aumentar el radio del electrodo, cabe 
A 
esperar que la reobase de voltaje presente un mínimo para un cierto radio del electrodo. Este 
comportamiento de la reobase de voltaje, sobre el cual Werner Irnich llamó la atención por pri- 


mera vez, se ha podido confirmar en trabajos experimentales. 


Si los electrodos son impolarizables, K (2) y L(s) se anulan al igual que pls). 


En ese caso las curvas amplitud umbral versus duración se comportan igual para pulsos de co- 
rriente y para pulsos de voltaje. Pero si los electrodos se polarizan, de la fórmula (18) se des- 
prende que cuando la duración del pulso aumenta, las amplitudes umbral de los pulsos rectangu- 
lares de voltaje se aproximan a su reobase más rápido que las amplitudes umbral de los pulsos 
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rectangulares de corriente se aproximan a su propia reobase. 
Por definición, las constantes de tiempo para los pulsos de voltaje y de corriente vienen dados 


lim, p(2, U, (1, )) sa _lim, ot, 1, (1,)) 


por los límites Ñ e A 
Rb Rb 


lsu 


Entonces de (5) y (18) resulta, teniendo en cuenta que pl0) = 7] : 
R 


Pp 


a ls _ R, 
so 1+p(0) R,+R, "sa pal 
Los resultados experimentales confirman las predicciones que se desprenden de la fórmula (30) 
para los electrodos convexos convencionales (Irnich, 1980): en igualdad de las demás condicio- 
nes, si los electrodos se polarizan, la constante de tiempo para los pulsos umbral de voltaje de 
forma rectangular es siempre menor que la constante de tiempo para los pulsos umbral de co- 
rriente de forma rectangular, y la diferencia entre ellas aumenta cuando aumenta el valor de la 


resistencia de la interfaz R, en condiciones de corriente continua. 
Además, para un intervalo de radios de electrodo tal que R p es un orden de magnitud mayor 
que R, , ambas constantes de tiempo deberían aumentar con el radio en forma análoga. 


La Figura 11 presenta un esquema idealizado de un circuito para estimulación monopolar. Está 
formado por el generador de pulsos del marcapaso, el cable del cátodo con una resistencia óh- 


mica KR. , un modelo de la capa doble electroquímica en la interfaz entre la superficie de la ca- 
beza del electrodo y los tejidos y la resistencia de los tejidos representada por dos resistencias 
óhmicas en serie, Re (correspondiente al casquete de tejidos no excitables) y R, 7 (correspon- 


diente al miocardio y el resto de los tejidos hasta alcanzar el ánodo ubicado en la superficie de la 
unidad generadora de los pulsos). 


Kc 
o 
CAPA DOBLE 


ELECTROQUIMICA 
FUENTE DE 
CORRIENTE 
2) REGULADA 


Figura 11. Esquema idealizado del circuito formado por el marcapaso y el área bioeléctrica pasiva. 


FUENTE DE 
POTENCIAU 
REGULADO 


] TEJIDOS 


BIOLÓGICOS 


En la figura se puede ver dos fuentes en el generador de pulsos: una de voltaje controlado y la 
otra de corriente controlada. La capa doble electroquímica está formada por una fuente de fem 
de corriente continua que representa la diferencia de potencial galvánico de equilibrio, a través 
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de la interfaz, en serie con la combinación en paralelo de un capacitor de Helmholtz C=C, 


con un resistor farádico Z=R;. 


En este caso particular: 


R,=R,+Ry¡ +R;, (31) R,=Z=R, (32) 
x(0=2exp Ll (20) 63) L(0)= exp 10120) 64) 
Es Tp Tu Ty 
e 
=REG 35 === 36 
Tr FH.H ( ) Ty Tr “E E: ( ) 
R, 1 
O EE 37 
ol) R 00 z (37) 


Para un modelo de cátodo esférico de radio 7, rodeado de un casquete interno de tejidos no 


excitables de espesor promedio e(n) (donde 1 representa una variable temporal discreta, por 


ejemplo los días consecutivos en los que se efectúan mediciones como las informadas en Fian- 


dra, Suárez-Antola y otros, 1981) y conductividad promedio O, (n) , y de un casquete externo 
infinito que comprende al resto de los tejidos, de radio interno 7; + e(n) y conductividad pro- 


medio O, , Se tiene, siendo (Y un factor de corrección geométrica: 


1 o 
R 5, (1) ATRAE (38 a) 


s 47.05, (n) A +2(n) 


1 1 


R zz ———| ——— 
ni) 470.0, 1, +2(n) (289) 


Las fórmulas (38) permiten describir la variación de KR, =R +R7,¡+R 7, asociada a las modi- 


ficaciones que sufren las propiedades de la barrera eléctrica interpuesta entre la superficie del 
cátodo y el tejido excitable, una vez que se produce la reacción inflamatoria al cuerpo extraño 
luego de implantado el electrodo. 

El edema, la infiltración celular y la formación de una cápsula fibrosa van alterando la conduc- 
tividad y el radio externo del casquete de tejidos no excitables adyacente al electrodo. 


Las variaciones en el espesor e(n) del casquete de tejidos no excitables y en R h (n) , Junto 


con las fórmulas halladas para las variables globales del umbral, suministran una base teórica 
para analizar el incremento fisiológico del umbral. 


La distancia promedio € aumenta progresivamente con el desarrollo de la reacción tisular, pero 
una vez completado y retraído el depósito de fibrina, € se estabiliza, generalmente después de 
haber sufrido una cierta disminución (Szabo y Solti, 1975). 

Si se produce una infección asociada al proceso reactivo a micro traumas mecánicos causados 
por la motilidad ventricular, las variaciones en la conductividad eléctrica y en la distancia entre 
el electrodo y el miocardio pueden ser aún más considerables. 
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Discusión y Conclusiones 


-Para estudiar los problemas que se plantean en relación con los umbrales de estimulación eléc- 
trica del miocardio parece conveniente introducir una distinción entre excitabilidad global y 
excitabilidad local en el marco de una descripción funcional del sistema electrodo-miocardio 
como la propuesta. 


- La generalización de un modelo de cable tridimensional excitable y no lineal, para describir la 
polarización de las membranas de un medio anisótropo y homogéneo, permitió dotar de un fun- 
damento teórico al concepto de masa crítica introducido por Lindemans para describir el umbral 
de estimulación del miocardio mediante electrodos externos. Se obtuvo una ecuación que res- 
tringe las dimensiones posibles de esa masa crítica, asimilada a un elipsoide (o a un hemi- 
elipsoide) crítico de células miocárdicas. 


-Se modeló el tejido cardíaco como un medio conductor homogéneo y anisótropo, compuesto 
por dos medios superpuestos (el intracelular y el extracelular) conectados por membranas exci- 
tables. No se tuvo en cuenta los posibles efectos de la anisotropía desigual ni de heterogeneida- 
des en los patrones de polarización producidos por el cátodo. El presente modelo del flujo de 
corriente en el miocardio ventricular, en su estado actual, parece demasiado simplificado. 


-En el modelo de membrana utilizado en este trabajo no se tuvieron en cuenta las variables de 
recuperación. Para continuar con el abordaje analítico, se podría incorporar una variable de re- 
cuperación mediante un modelo matemático del tipo del introducido por Fitz-Hugh y Nagumo 
(Fitz-Hugh, 1969). 


-A diferencia de la longitud límite de Rushton deducida por Noble para una fibra excitable, la 
deducción de las dimensiones de los elipsoides críticos efectuada en este trabajo muestra que a 
esos elipsoides les corresponde una familia de infinitos estados de polarización umbral que pue- 
den aparecer en un mismo miocardio. 


-Cual sea la región crítica que en cada caso corresponda tener en cuenta para estudiar la génesis 
de un potencial de acción, depende de las regiones de influencia de los electrodos sobre el mio- 
cardio. Esas regiones se caracterizan por un apartamiento significativo de los voltajes trans- 
membrana de sus valores en reposo, provocado por el campo de corrientes eléctricas generado 
por el electrodo del marcapaso. 


-Para estimulación mediante pulsos rectangulares de corriente controlada, las fórmulas obteni- 
das predicen el aumento de las reobases de corriente y las constantes de tiempo para pulsos de 
corriente, con el incremento del radio del electrodo, el incremento del espesor del casquete de 
tejidos no excitables interpuesto entre el electrodo y el miocardio y con el cambio en la modali- 
dad de estimulación de monopolar a bipolar. 


-Para estimulación mediante pulsos rectangulares de voltaje controlado, las fórmulas obtenidas 
utilizando un modelo lineal para describir las propiedades eléctricas de la interfaz electrodo- 
tejidos, relacionan las reobases de voltaje y las constantes de tiempo (para pulsos de voltaje 
controlado) con la resistencia variable de los tejidos, la resistencia de la capa doble electroquí- 
mica de la interfaz a la corriente continua, y con las reobases de corriente eléctrica y las cons- 
tantes de tiempo para pulsos de corriente, respectivamente. 

Predicen las modificaciones en el umbral de voltaje con las variaciones en las dimensiones, 
forma y polarización del electrodo, así como con las modificaciones en la resistencia de los 
tejidos asociadas a la formación del casquete de tejidos no excitables. 
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-Empleando un modelo de cátodo esférico con un factor geométrico de corrección para tener en 
cuenta geometrías más realistas, y asumiendo un comportamiento óhmico para el conductor de 
volumen de los tejidos en las condiciones imperantes durante la estimulación del corazón me- 
diante electrodos externos, se obtuvieron fórmulas que relacionan las resistencias del casquete 
de tejidos no excitables y la del resto de los tejidos, con las conductividades eléctricas promedio 
y el espesor promedio del casquete de tejidos no excitables. Esas fórmulas permiten describir 
los efectos de la reacción al cuerpo extraño que se produce una vez implantado el electrodo. 


-Los eléctricos modelos lineales para las interfaces entre los electrodos y los tejidos pueden 
resultar en una simplificación excesiva de la complejidad real de la interfaz. El umbral de no 
linealidad depende, entre otras variables posibles, del material del electrodo y su forma (Jaron, 
Briller y otros, 1969). 

Una vez franqueado ese umbral, se modifica el capacitor de Helmholtz (generalmente aumenta) 
al igual que la fuerza electromotriz que surge de la interacción entre varios procesos fisicoquí- 
micos que se desarrollan en diferentes lugares y en diferentes instantes. La resistencia farádica, 
considerada como parámetro constante, deja de ser aplicable. 


- Si bien la evidencia experimental conocida por el autor confirma las predicciones que resultan 
de las fórmulas obtenidas en este trabajo, en general esa confirmación es solamente cualitativa. 
Esto es debido a que la evidencia proviene de experimentos que no se diseñaron específicamen- 
te para poner a prueba las fórmulas presentadas y discutidas en este trabajo. 


-Queda entonces para investigaciones posteriores deducir y contrastar las consecuencias que se 
obtienen del abordaje teórico presentado aquí, corregir los errores cometidos y sustituir el pre- 
sente marco teórico por uno más adecuado para el análisis de la estimulación del corazón me- 
diante electrodos de marcapaso. 


Apéndice: 


Consideremos ahora, brevemente, el proceso que conduce a la fórmula que permite determinar 


si un pulso de corriente catódica de duración 1, es un pulso L, (e ) justo umbral: 
5 . 
f: G, (1, -1)i, (1) dt=0, 


Se puede partir de un replanteo que generaliza una observación debida a Trautwein (Trautwein, 
1975): para que la estimulación sea exitosa es necesario transportar una cantidad crítica de carga 
a través de un área adecuada de membrana celular excitable, membrana que se encuentra ini- 
cialmente en reposo, tal que permita que esa membrana se polarice hasta un estado umbral. 
Dicho de otro modo, el electrodo estimulador debe ser capaz de generar en el lugar adecuado un 
patrón de densidad de corrientes locales de suficiente amplitud y duración como para que se 
produzca una transferencia de carga a través de las membranas que supere un valor crítico (Foz- 
zard, 1979). 


Suponemos que el cable tridimensional ocupa la región z > 0. Introduciendo las constantes de 


E R R : 
espacio 22 =“% 22 =% 22 =“2 y la constante de tiempo 7, =R C, .,, la ecuación 
xx R y y R ZZ R m m mV 
Xx y z 


del cable tridimensional (9) linealizada en torno al voltaje de reposo se puede poner de la forma 
(refiriendo los voltajes al de reposo): 
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Asumiendo ahora que, para z constante, los perfiles de despolarización en las nuevas variables 
espaciales verifican V,, (1 X.Y, z) =V, ( AY, ¿) , la Figura 6 (bis) siguiente muestra, 


para un instante de tiempo dado, dos perfiles de potencial transmembrana correspondientes a 
dos valores de z. Ambos presentan simetría de revolución respecto del eje z . Cerca de éste son 
cóncavos y más lejos son convexos. 


| 


Figura 6 (bis). La figura muestra la despolarización de las membranas del modelo de cable tridi- 
mensional situadas en dos planos ortogonales al eje de simetría. Uno más próximo y el otro más 
alejado respecto del cátodo. 


La Figura 12 muestra, a la izquierda, la traza de uno de esos perfiles de voltaje transmembrana 
en un plano que pasa por el eje de simetría. Por la simetría de revolución asumida, esas trazas 
son todas iguales al pasar de un plano por el eje a otro, siempre que z se mantenga constante. 
Aproximamos la superficie por su paraboloide osculante en el punto X =Y =0, como aparece 
indicado a la izquierda. 

A la derecha aparece la traza circular de ese paraboloide en el plano de las variables X, Y 
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Van 


A 


Erazo ele? 
perfil e, 
f LS 


des pol: jár Jue dt 


A 1, , 
£rasa drA y 
paróbo fa ¡ele 1 


oscufonte 
J 


| 
5 E 


Dep = >. - 
1O ¡8 PS 

PA que se 

emplts en a 


en elef 


ro dp reo mas 
QA Ñ yA dues cn. ", 
py ed, le anvta € 


parabo/riee osevlonte) 
Figura 12. Para un mismo valor de 7 : A la izquierda se esbozó la traza del perfil de voltaje trans- 
membrana en el cable tridimensional y la traza del paraboloide osculante a ese perfil. A la derecha, 
se observa la traza circular del paraboloide osculante sobre el plano determinado por las variables 


> AO, CE 
A As, 
Como los voltajes se refieren ahora al voltaje de reposo V,, es posible aproximar 
o o As, 
5 V,, (1,0,0,2) + 5V,, (1,0,0,2) por -y 5 V, (+,0,0, 2). 
OX OY (1, +eY 


Esta aproximación se obtiene comenzando con el desarrollo de Taylor del perfil de voltaje, has- 
ta el segundo orden en torno al punto X =Y =0, para t y z dados: 


V, (1,X,Y,2)=V,(1,0,0,2)+X — (1,0,0,2)+ y As (1,0,0,2)+ 


m 


(1,0,0,2)+0(X?+Y?) 


Ly Velro0z3xr a Va —— (1,0,0, z)+ ¡Ly 
A). OX OY DOY 


Si se desprecian los términos de orden superior al segundo se obtiene la ecuación del parabolide 
osculador. El término que se presenta en las derivadas cruzadas se anula por razones de sime- 
tría. 

Si P'representa el círculo donde este paraboloide corta el plano X, Y, integramos sobre ese 
círculo y obtenemos: 


Jv, (2, X,Y,2) ds = (1) (1,0,0,2)+2 as Pa =.(1,0,0,2) +5 ES ds a =.(1,0,0,2) 
Jr i 212 ox? Y 


Como V,, (z, X,Y, z) se anula sobre I”, el miembro de la izquierda es nulo. 


Además: fx0as=[Y?as=2 JL0+Y)as 
E: r 2 T 


1 
Si R es el radio del círculo I”, ((1)=27R y > [10 +7?) ds = R' por lo cual: 


Tr 


2 2 
(1,0,0,2)+ =— (1,0,0, 2) = V,, (1,0,0, 2) 


47 o OA E O 
Asumimos que, independientemente de z, =>=Y ———5G siendo ——= + 
O A O 


Entonces, sobre el eje de simetría: 
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2 
—V (1,0,0 Z)+ 1+y == 
m m A: 2 m 2 Z a 2 
Ot (r, + e) 550 Z 
Integrando entre z=0 (el punto más próximo al cátodo) y Z=>+00 , teniendo en cuenta que 


0 
—V (2,0, 0, z) debe tender a cero cuando z—>-+00, y multiplicando ambos miembros por la 
L 


capacidad C,, y por unidad de volumen del cable tridimensional, definiendo la carga por unidad 


de área (1) = Cry [val a £,0, 0,z) d z, la constante de tiempo f,, =———">5— y la 


A, Am lrmte) +e) 


dra +(m+ey +(m te) 


la nueva variable de estado o(t ) que resume la dinámica de la polarización): 


d 1 A O 
— — ez |-C ,—V (1,0,0,0 
PS el Lo js OZ a ) 


constante de espacio 4, = resulta esta ecuación diferencial ordinaria (para 


1) ot a 
Por otra parte: ——V,, (2,0, 0,0) e E donde O es la conductividad del con- 
Oz 4ra(r+e) O 


ductor de volumen de los tejidos. 
Integrando la ecuación para o(t ) entre 1=0 y t=1 p “on (0) =0 (pues trabajamos llaman- 
do V,, alo que denominamos V,, —V, en el cuerpo principal del texto) se obtiene: 


[exp le, -+) i(e)di=C, dra 0 (1 +e) o(t,)=0(1,) 


0 Lor 2Z 
Si el pulso catódico es un pulso justo umbral 7, 0) cuando al final del pulso la carga por unidad 


de área o(t, ,) alcanza un valor crítico (, , se obtiene la condición: 


les 0-1) za) 


0 Lol 


¡, (1) dt =0, (A1) 


Donde por definición aparece una cantidad de carga crítica: 


T 
0. ==" 4Ara 0 (r,+e) o. (A2) 
aj A. , ) 
Entonces: 
O  4rac ¿ 2 
Li = | A Hr +e o, (A 3) 
Rb he - Cu Az Y y ( 0 ) s 
Finalmente: 
To 
L.¡= AA (A4) 
a 
(7, +e) 
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Con esto se justifica el marco teórico planteado en las páginas 12 y 13, en la sección “Deduc- 
ción preliminar de la fórmula de Lapicque-Hill: análisis de los umbrales de captura mediante 
pulsos de corriente controlada.” 
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